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STABILITY DERIVATIVES AT SUPERSONIC SPEEDS 
OF THIN RECTANGULAR WINGS WITH DIAGONALS AHEAD OF TIP MACH LINES 

By Sidney M. Hahmon 


SUMMARY 

Theoretical results are obtained, by means of the linearized 
theory, for the suif ace-velocity-potential functions, surface- 
pressure distributions, and stability derivatives for various 
motions at supersonic speeds of thin flat rectangular wings 
without dihedral. The investigation includes steady and accel- 
erated vertical and longitudinal motions and steady rolling, 
yawing, sideslipping, and pitching for Mach numbers and 
aspect ratios greater than those for which the Mach line from 
the leading edge of the tip section intersects the trailing edge 
of the opposite tip section. The stability derivatives are derived 
•with respect to principal body axes and then transformed to a 
system of stability axes. In the case of yawing, a treatment for 
the infinitely long wing which takes account of the spanwise 
variation in the stream Mach number is extended to the finite 
wing, and a plausible, although not rigorous, solution is obtained 
for the wing tip effects. 

The results for this investigation showed that positive yawing 
at supersonic speeds may produce a negative rolling moment 
•in contrast to the behavior at subsonic speeds where a positive 
rolling moment would be produced. The attainment of super- 
sonic speed should produce a significant change in the positive 
direction of the yawing moment per unit rolling velocity. The 
results also indicate that unstable tendencies are produced by 
vertical accelerations if 

where A is wing aspect ratio and M is stream Mach number. 

INTRODUCTION 

Recent developments in supersonic airfoil theory (refer- 
ences 1 to 4) have led to the calculation of many of the 
supersonic stability derivatives for various plan forms. In 
references 5 to 8, various theoretical supersonic stability 
derivatives for small disturbances are presented for thin flat 
wings of delta plan form. In reference 9, the supersonic 
damping due to rolling is given for triangular, trapezoidal, 
and related plan forms. 

In the present paper the methods of references 4, 10, and 
11, which are based on the linearized theory for a unif orm 
stream Mach number, are used to derive the supersonic 
surface-velocity-potential functions for thin flat rectangular 
wings without dihedral in steady and accelerated vertical 
motions and steady rolling, sideslipping, and pitching 
motions. The potential functions that are obtained are then 


used to derive formulas for the pressure distributions and the 
stability derivatives for the foregoing motions and also for 
steady yawing. In the case of yawing, a simple treatment 
given in reference 7 for the infinitely long wing, which takes 
account of the spanwise variation in stream Mach number 
associated with yawing, has shown that the assumption of a 
uniform Mach number is far from adequate to describe the 
compressibility effects. This treatment is extended herein 
in order to evaluate the wing tip effects for the yawing finite- 
span wing. 

The steady motions that are treated herein are assumed 
to give small deviations from the undisturbed flight path 
and the accelerated motions are assumed to have small 
accelerations. Theoretical results based on this assumption 
for steady motions have, in general, been found to be reliable; 
however, the reliability of such results for unsteady motions 
is as yet unverified. The results presented herein cover a 
range of Mach number and aspect ratio greater than that 
for which the Mach line from the leading edge of the tip 
section intersects the trailing edge of the opposite tip section. 
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rectangular coordinates (see fig. I) 
induced flow velocities along x- and y-axes, 
respectively 

coordinate in flight direction if this direction 
is inclined to x-axis 

incremental flight velocities along x-, y~, and 
z-axes, respectively (see fig. 2) 
derivative of u with respect to time 
accelerated vertical motion 
undisturbed flight velocity 
local flight velocity after disturbance; used to 
indicate inclination of flight direction to 
x-axis (see fig. 1) 

angular velocities about x-, y-, and z-axes, 
respectively (see fig. 2) 
speed of sound 
stream Mach number {Vja) 

Mach angle ^sin -1 -^^ 

wing angle of attack in steady flight, radians 

(w/T) _ 

local inclination of airfoil surface with respect 
to free stream, radians 
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C x 
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L,M,N 


derivative of a with respect to time 
time following disturbance, seconds 
angle of sideslip, radians (v/V) 
chord 

wing semispan 
wing span 
total wing area 

region of integration over portion of wing 
sin-face (see fig. 3) 


aspect ratio 



distance of origin of stability axes from the 
midchord point, measured along 2 -axis, 
positive ahead of midchord point 

mass density of air 

disturbance-velocity potential on upper sur- 
face of airfoil 

auxiliary variables which replace x and y, re- 
spectively (see fig. 1) 

indicates a transformation of origin of x- and 
y-axes or £- and ij-axes from leading edge of 
center section to leading edge of tip section 
( ya—y—h on right half-wing; y a =—y—h 
on left half-wing) 

pressure difference between lower and upper 
surfaces of airfoil, positive in direction of 


lift 


nondimensional coefficient expressing ratio of 
pressure difference between lower and 
upper surfaces of airfoil to free-stream 


dynamic pressure 




constant given by equation (9) 
induced suction force on wing tip per unit 
length of tip 

forces parallel to z-, y~, and 2 -axes, respectively 
(see fig. 2) 

/ X 

longitudinal-force coefficient ( - 


lateral-force coefficient / - 




V 2 SJ 


\§ v 2 s 

(z 

— 

f ™ 

skin-friction drag coefficient 
'"Skin-friction drag 


) 


C 


§ V 2 S 


') 


moments about z-, y-, and s-axes, respec- 
tively (see fig. 2) ; M is also used to refer to 
Mach number 


C, 

c m 

c n 

Subscript: 

1,2 


Superscript: 


rolling moment coefficient 



pitching-moment coefficient 



yawing-moment coefficient 



contributions of normal pressures and skin 
friction, respectively, to C Hf ; also used to 
indicate component parts of ft* c a&t 
Cz u> and 


w,p contributions caused by vertical motion and 

rolling motion, respectively 

Whenever p, g, r, /3, «, u, a, and u are used as subscripts, 
a nondimensional derivative is indicated and this derivative 
is the slope through zero. For example: 




Unprimed stability derivatives refer to pi-incipal body 
axes; primed stability derivatives refer to stability axes. 


ANALYSIS 

GENERAL CONCEPTS 

The coordinate axes and the symbols used in the analysis 
of the rectangular wing are shown in figuro 1. The derivation 
of the formulas for the surface-veloeity-potential functions, 
pressure distributions, and stability derivatives is made 
initially with reference to principal body axes which aro fixed 
in the wing with the origin at the midchord of the center 

section 0, 0^. This system of axes is shown in figure 2 (a). 

The transformation of these stability derivatives to a system 
of stability axes (fig. 2 (b)) is discussed in the section entitled 
"Results and Discussion.” 

The stability derivatives are determined from integrations 
of the forces and moments over the wing. For vertical and 
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pitching motions which yield equal and opposite suction 
forces along the edge of each wing tip, the only resultant 
forces and moments acting on the wing, if skin friction is 
neglected, are those caused by the pressures on the airfoil 
surfaces. These pressures are obtained from the familiar 
Bernoulli equation. In rolling, yawing, or sideslipping, how- 
ever, unbalanced suction forces which produce lateral forces 
and yawing moments are induced along the wing tips in 
addition to the forces and moments resulting from the pres- 
sure normal to the wing surface. The subsequent analysis 
for the calculation of the stability derivatives is then resolved 
to a determination of the pressure distribution normal to the 
surface and the unbalanced suction forces along the wing-tip 
edges. 

The pressure difference between the upper and lower 
surfaces (positive upward) at any point on the wing is 
determined from the general Bernoulli equation for small 
disturbances as 

iP-2 P (r'|f+!) (i) 



(a) Principal body axes. Origin at center of wing. 

(b) Stability axes. Origin at point x«,,0,o). 

Principal body axes dotted tor comparison. 

FiGCEZ 2. — Velocities, forces, and moments relative to principal body and stability axes. 

where V is the local flight velocity and l refers to a coordi- 
nate measured in the flight direction. The term 5<£/5f ex- 
presses the effect of any unsteadiness in the flaw. The 
velocity potential <f> in equation (1) must be determined so as 
to satisfy the linearized partial-differential equation (with 
time dependency if the motion is unsteady) of the flow and 
the boundary conditions associated with the particular mo- 
tion under consideration. Thus, the potential must give 
streamlines that are tangent to the airfoil surface and a pres- 
sure field that is continuous at all points exterior to the 
wing. Equation (1) shows that the pressure distribution 
on the wing is determined when the surface-potential function 
is found. 

The method of referenced is in general adaptable to the prob- 
lem of obtaining the surface-potential function <p in supersonic 
flight to meet boundary conditions associated with small 
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steady motions, such as vertical motions, rolling, sideslipping, 
and pitching. The method is an extension, to include tip ef- 
fects, of the work of Puckett and others which uses the super- 
position of elementary source solutions to obtain the potential 
function. In cases where a point on the wing is influenced 
by two or more mutually interacting external fields, the inter- 
action introduces difficulties in the solution for the surface 
potential. (See also reference 12.) If any point on the wing 
is influenced by only one independent external field, however, 
the potential function in a region affected by the wing tip 
may be obtained by integration of elementary source solu- 
tions solely over an appropriate area of the wing. The 
strength of these sources is shown to be a function only of the 
local slope of the airfoil surface with reference to the free- 
stream direction. Inasmuch as the slope of the airfoil sur- 
face with reference to the free-stream direction is known for a 
given motion, the distribution of sources is known and, con- 
sequently, the distribution of the surface-potential function 
is determined by an integration of the elementary source 
solutions over an appropriate area of the wing. 

As applied to the rectangular wing at supersonic speeds, 
the foregoing method of reference 4 for one independent 
external field is valid as long as the foremost Mach wave 
from one tip does not intersect the opposite tip, that is, for 
Mach numbers and aspect ratios for which ABsS 1. For 
this case, the potential at a point on the top surface of a 
tlfin flat wing may be determined by means of equation (14) 
of reference 4 and is as follows: 




d£ dr] a 


^W-B^Va-Va ) 2 


( 2 ) 


where a' represents the local angle of attack of the airfoil 
surface at the point (£,ij„). Figure 3 shows a typical region 
S w for determining the potential at a point (x,y a ) in a 
rectangular wing. The figure shows the boundaries S a over 
which the integration must he performed, for a point (x,y a ) 
which is affected by the wing-tip region. If the point 
(x,y a ) is located at or inboard of the foremost Mach line 
from the tip, this point is unaffected by the tip region and 
S w is bounded by the leading edge and the Mach forecone 
from (x,y a ). Suppose that the surface potential <t>(x,y ) has 
been obtained from equation (2) or by some other method, 
then the differentiation of $ with respect to the coordinate 
in the free-stream direction determines the pressure distri- 
bution by means of the Bernoulli relation, equation (1). 

The expressions for determining the surface potential and 
the pressure coefficient for unsteady motions are discussed 
in the section entitled “Derivation of Formulas." 


DERIVATION OF FORMULAS 


along the wing tips. The integrals required for these deriva- 
tions and also those required for the stability derivatives am 
integrablo either directly or after reduction by parts by means 
of standard formulas such as are given in reference 13; 
hence, the details for the integrations are not shown. 
In the operations involving factoring from radicals, care 
must be used to preserve tho correct sign of the factors; for 
example, if 

y*< o 


then 


v^?=V(-y«)*=-y« 


For brevity, the final formulas are omitted from the deriva- 
tions and appear only in tables at tho end of the paper. 
Thus, the distributions of <f> and A c P are summarized in 
table I, and the stability derivatives are summarized in 
table II. 

All the derivations arc made specifically for a whig for 
which AB$z 2, that is, for winch the foremost Mach wave 
from a tip does not intersect tho remote half-wing. The 
formulas in table I for the potential <f> and pressure coefficient 
Ac P that are obtained for AB s^2 can be applied to wings in 
which 1 ^ AB ^ 2 by using tho principle of symmetry and 
superposing separately each tip effect at tho point under 
consideration to the value obtained for tho infinitely long 
wring. A consideration of this superposition principle for the 
rectangular wing shows, however, that tho stability deriva- 
tives which are obtained for AB>. 2 apply as well to wings 
for which AB^ 1. A more detailed description of table II 
is given in the section entitled “Results and Discussion.” 


VERTICAL. PITCHING. AND LONGITUDINAL MOTIONS 


Derivatives —Cz a and — C Zj . — For steady pitching motion 
about a lateral axis tlirough the midchord point, the local 
slope of the airfoil surface with respect to the free-stream 
direction is 


a’— a- f- 



where a is the angle of attack in the absence of pitching. 
In order to obtain the potential distribution, this value of 
a' is substituted into equation (2) and tho doublo integration 
for the variables £ and i ?„ is performed between the limits 
indicated in figure 3. The pressure coefficient is then ob- 
tained from equation (1) for steady motions as 


Ac P 


4 4 5<£ 

-yUo-Vbx 


( 3 ) 


The subsequent derivation of formulas for the various 
motions will involve first the determination of distributions 
of surface potential and then the determination of surface- 
pressure distributions and any unbalanced suction forces 


These pressure coefficients are then differentiated with respect 
to a and g. The integrations of the respective distributions 
of Ac P over the wing and conversion to nondimensional 
units then give the derivatives — C z<x and —C Zf . 
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Wing tip 



FlOl'BI J.— Region of integration Cor obtaining supersonic velocity-potential function for 
rectangular wing of Unite wing span. Region of integration Cor point shown cross- 
batched: r«=g-* on right half-wing; i r—k on left half-wing. 


Derivatives C m and C- . — Derivatives C m and C m 
are obtained directly from the pressure-coefficient distribu- 
tions for angle of attack and pitching, thus 

where A c P for angle of attack and pitching is given in table I. 

Derivatives C Xa and C Xf . — At supersonic speeds the re- 
sultant pressure force on a rectangular wing of zero thickness 
acts normal to the surface as there is no suction at the 
leading edge. Thus, the forces in the x-direction arise solely 
from skin friction. On the assumption that the skin friction 
is independent of a and q, the derivatives C Za and C Xt are 
zero. 

Derivatives —C Zu , C„ u , and —C Zu . — The derivative —C Zu 
is obtained from the equation 

-z=-a r c Za % (v+uys 

Then 

~ Cz =V^Z l-a'Cz.ty+u)^ 


The function— C Za is obtained from table II. Its depend- 
ence on the incremental flight velocity U is indicated by 


g B in the form — 


I and 


W 

T+if’ 


Then 



The derivative C„ a is obtained from the equation 


Then 


M=a'C Ua ?(V+u)*Sc 
C*u=y [a / C 0 , a ('r+'u) s ] B -* ) 


The function is obtained from table II, where its de- 
pendence on u is indicated by writing a' and B in the same 
form given previously for the derivative —C a ; u - Then 


w(V-\-u) 




The derivative ~Gx u results from skin friction and is 
obtained from the equation 


-X=C Dq f (F+u)*S cos « 

Then 

ACCELERATED MOTIONS 

For accelerated motions in the vertical plane of symmetry, 
the pressure coefficient from equation (1) is 



The surface potential <j> in equation (4) for unsteady motions 
of thin airfoils in two-dimensional supersonic flow has been 
derived in reference 14. In reference 10, the methods of 
references 4 and 12 for steady flow at supersonic speeds are 
extended in order to determine solutions for the surface po- 
tential and pressure coefficient for unsteady motions in three- 
dimensional flow. In the present analysis the solutions 
obtained in reference 10 will be utilized to calculate the deriv- 
atives in vertical motions with small constant accelerations. 

Derivatives —C Zi , C ma , and —C z& . — The surface potential 
4> for uniformly accelerated motion as obtained for the region 
within the tip Mach cone is (reference 10, equation (31)) 

(ra*-^^cos-'(^+l)} (5) 

In converting from the notation of reference 10 for a 
rectangular wing to the present notation, the following 

V 

transformations are made: U=V, m=a, a=0, (3=5, c=^-> 
h=l, Ui=v u =^ (x+y s B), (x—y a B),^= 1, and 
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C v = — In order to obtain in the region between the tip 

Mach cones, y a in equation (5) is set equal to — -g and, 
therefore, 

, V { . x APax > \ 

4>— g 2W / ^ 

The pressure coefficient Ac P contributed by the vertical 
accelerating motion is obtained by partial differentiation of 
4> (equations (5) and (6)) with respect to x and t, by letting 
t=0, and then by substituting these expressions for b<t>/bx and 
b(j>fdt in equation (4). This process yields in the region 
within the tip Mach cones __ . ... 

Ac P (x,y a ) ("^ 1 )+ 2 -S 2 ^/-y a (y«+g)] 

■ (7) 

and in the region between the tip Mack cones 

A C P = (8) 


from suction along the tips. These suction forces may be 
evaluated by applying a method suggested in reference 15 
for incompressible flow and modified for compressibility 
effects in reference 3. Thus, if the induced surface velocity 
normal to the wing tip is expressed as 

Po=±-r== I (9) 

V — 2/o 

where Q is a constant, then the suction force per unit length 
of tip is 

F S =TTp(P 

(A more general expression for edge suction dial is still valid 
when the edge is inclined to the stream is given in reference 3 
and recast in reference 7.) 

Consider the induced surface velocity normal to the tip of a 
wing rolling with an initial angle of attack a. This velocity 
is 


Equations (7) and (8) correspond to equation (33) of 
reference 10, after the appropriate transformations noted 
previously for <f> are made. 

The derivatives -G z& and C m& are then obtained by 
integration of the corresponding Ac^-distributions and con- 
version to nondimensional units. The derivative Gx & is 
shown to be zero by the use ol assumptions similar to those 
noted previously for 0 Xa . 

Derivatives -C Zi , C mi , and -Cxy — For small accelerations 
along the flight path, the potential will remain substantially 
unchanged. The increments in pressure caused by these 
accelerations, therefore, are assumed to be negligible, and the 
derivatives ~ C H> C m<t , and -C x& are approximately zero. 

BOLLING 

Derivative Ci p . — In steady rolling motion with angular 
velocity p, the local slope of the airfoil surface with respect 
to the flow direction is 


where the superscripts w and p refer to the potentials obtained 
for a vertical motion and a rolling motion, respectively. 
From table I 

<i> u (x,ya) [§ cos- 1 (^+ 1 )+ 2 -\/- 2 /«(y«+|)] 
and partial differentiation of <j>° with respect to y yields 



where 0. Very near the tip, -y a -*0 and 



The potential in rolling <j> p is given in table I. By partial 
differentiation of 4> p with respect to y and then by letting 
-y a —» 0, there results 


pv jpbh±]>) 
V~ V 


In order to obtain the potential distribution this value of a' 
is substituted into equation (2) and the double integration 
for the variables £ and is performed between the l i mits 
indicated in figure 3. The pressure coefficient is then ob- 
tained from equation (3). The derivative C lp is obtained 
by integrating the moments of the Bernoulli pressure dis- 
tribution for rolling given in table I and by converting this 
result to coefficient form. . . 

Derivatives C Yp and C np . — In a rolling motion, the lateral 
force and yawing moment relative to body axes result entirely 



The resultant induced surface velocity normal to tho wing 
tip as —y a ~* 0 is obtained by adding equations (11) and (12). 
Thus 



(13) 


The plus sign before V refers to the right wing tip and the 
negative sign refers to the left wing tip. 
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Very near the wing tip, equation (13) has the same form 
as equation (9) and, therefore, the total suction force per 
unit length along the wing tip is 


[y (»-£,)+ w] 

(14) 

In equation (14) only the term ± ^ pX ^ a will 

give rise to a lateral force and a yawing moment which are 
obtained by integrating this term along the wing tips. 
These forces and moments are then converted to non- 
dimensional form to give the derivatives C Tp and C Xp . 

SIDESLIP 


The pressure coefficient obtained from equation (1) for 
steady flight is 


4 b<j> 
ACp ~Y' bl 


where V' and l are measured in the flight direction. If 
sideslip occurs the flight direction is inclined relative to the 
r-axis by the sideslip angle fi. The rectangular wing in 
sideslip, therefore, becomes equivalent to a yawed wing 
with the leading wing tip raked out and the trailing wing tip 
raked in. If the Kutta-Joukowski condition at the trailing 
wing tip is neglected, the potential function for the yawed 
rectangular plan form may be obtained by the method of 
reference 4. In reference 11, the method of reference 4 is 
extended in order to obtain solutions for edges for which the 
Kutta-Joukowski requirement must be satisfied. 

Physical considerations suggest, however, that for small 
sideslip angles, the actual flow for typically rounded wing tips 
would in general be unlikely to conform to the Kutta- 
Joukowski conditions along the trailing wing tip. The edge 
suction for a lifting wing arises because of the flow from the 
bottom surface to the top surface around the side edge. This 
flow may be presumed to go around any boundary layer that 
may be present. The local boundary layer thus experiences 
the edge suction. Rough calculations suggest that the edge 
suction per unit area is approximately constant from the 
leading edge to the point of maximum profile thickness, and 
then increases rapidly from the point of maximum thickness 
to the trailing edge. The pressure gradient is therefore 
favorable and the flow at the side edge is not expected to 
separate. This condition should persist for small or moderate 
amounts of sideslip until the additional pressure increment 
caused by sideslip produces a strong adverse pressure 
gradient. Further theoretical and experimental investigation 
is required to obtain quantitative results regarding these 
phenomena. On the basis of the foregoing considerations, it 
will be assumed in the present analysis that the Kutta- 
Joukowski condition is not satisfied along the trailing wing 
tip. The effect of satisfying the Kutta-Joukowski condition 
along the trailing wing tip in sideslip is discussed in this 
analysis and also in the section entitled “Results and 
Discussion.” 


Derivative C tfj . — The potential corresponding to a thin 
rectangular wing at an angle of attack and a finite angle of 
sideslip may be obtained from reference 4, equation (20). 
The corresponding pressure distribution may be obtained 
from reference 11, appendix C, equation (C4). These solu- 
tions from references 4 and 11 were simplified to the approxi- 
mate form for small angles of sideslip (/3*-d) and converted 
to the present notation with respect to axes shown in figure 1 . 
The distributions for <f> and Ac P caused by combined vertical 
motion and sideslip are given in table I. The regions for 
which these expressions for <f> and A c P are applicable are 
bounded by Mach lines with respect to the stream velocity 
V which is inclined to the r-axis by the sideslip angle (3. As 
noted previously, these expressions do not satisfy the Kutta- 
Joukowski condition along the trailing wing tip. As indicated 
in reference 11, however, the Kutta-Joukowski condition 
along the trailing wing tip merely cancels the radical term 
in the expression for A c P within the Mach cone from the 
trailing wing tip. 

A consideration of the foregoing ACp-distributions indicates 
that as a result of sideslip the lift within the Mach cone from 
the leading wing tip is increased, whereas the lift within the 
Mach cone from the trailing wing tip is decreased. A rolling 
moment is thereby produced. Furthermore, as a result of 
sideslip, the Mach lines are shifted toward the trailing wing 
tip, and this shift contributes an additional rolling moment. 
The magnitude of the rolling moment caused by sideslip is 
given in table II in terms of the nondimensional derivative 

Derivatives C Yp and C np . — The derivatives C Yp and C Hp 
can result solely from suction forces which are induced 
at the wing tips. These suction forces for sideslipping motion 
were evaluated by a method similar to that described pre- 
viously for obtaining C Tp and C Uf . The treatment for side- 
slip was based on the conclusion, noted previously, that the 
Kutta-Joukowski condition is unlikely to be satisfied for 
typically rounded wing tips at small angles of sideslip. The 
potential <f> for determining the induced velocity normal to 
the wing tip was obtained from table I. The resultant lateral 
force and yawing moment are given in nondimensional form 
in table II. 

TAWING 

In yawing flight, the stream velocity varies linearly along 
the span. This effect introduces variations of both dynamic 
pressure and compressibility effects along the wing span. 
The surface potential as expressed in equation (2) satisfies 
the linearized potential equation for a uniform stream Mach 
number, but is inadequate to account for the compressi- 
bility effects associated with a spanwise variation of stream 
Mach number. (See reference 7.) The case of the trapezoidal 
wing with tips cut off along the Mach lines (raked tips) 
was analyzed in reference 7. It was shown that the pressure 
distribution could be obtained by application of the Ackeret 
two-dimensional theory modified by using the local Mach 
number at each spanwise station as affected by the yawing. 
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Inclusion of the spanwise variation in Mach number was 
demonstrated to have a profound effect on the pressure 
distribution. 

The addition of suitable triangular tips to the aforemen- 
tioned trapezoidal wing converts it into a rectangular wing. 
The added tips lie wholly within the tip Mach cones and 
thus their addition does not alter the pressures on the trap- 
ezoidal portions. A rigorous solution for the pressures on 
the tip portions cannot yet be demonstrated. However, an 
expression that appears plausible has been obtained. This 
pressure distribution for the tip portions is derived by super- 
posing on the Ackeret pressure distribution, as modified by 
local Mach number, an appropriate function which fulfills 
the boundary condition for no pressure discontinuities in 
the region exterior to the wing. This function thus repre- 
sents the effect of the wing cut-off and is designated herein 
as the tip effect. The pressure difference A P at any point 
according to the Ackeret theory based on local Mach number 
is (reference 7) : _ 



Equation (15) shows that the pressure distribution for an 
infinitely long wing which has a steady yawing velocity r and 
vertical velocity w is expressed by two components. One of 
these components is proportional to w, is constant, and gives 
the pressure distribution contributed by an angle of attack 
in straight flight. The other component is proportional to 
wr, gives a linear antisymmetrical distribution with respect 
to y, and expresses the pressure distribution contributed by r 
yawing. 

It will be recalled that the solution for steady rolling, 
treated in a preceding section, resulted likewise in a pressure 
distribution proportional to y in the region between the tip 
Mach cones. The pressure distributions contributed by roll- 
ing and by yawing are thus proportional in the region be- 
tween the tip Mach cones. The wing cut-off is effected by 
canceling the disturbance pressures outboard of the desired 
tip location by means of a function that satisfies the bound- 
ary conditions on the wing. Because the two pressure dis- 
tributions to be canceled correspond in the yawing and roll- 
ing cases, the incremental pressure function or tip effect for 
each case evidently must reduce to forms which will have 
the same factor of proportionality in the entire plane of the 
wing outboard of the tip. It seems reasonable to assume, 
therefore, that for small yawing motions the two pressure 
distributions will also have very nearly the same factor of 
proportionality within the tip Mach cones. 

The proportionality constant between the pressure distri- 
butions for rolling and yawing motions may be determined 
by a comparison of the cases of rolling and yawing in columD 4 
of table I. The pressure coefficient per unit yawing veloc- 
ity is seen to be ajB 2 times the pressure coefficient per unit 
rolling velocity, or 

olV 

(ACp) yaaint B 2 'p (ACp)r*m»* (16) 

where equation (16) will apply over the whole wing. 


Derivative C (r . — The preceding analysis indicated that the 
pressure -distribution per unit yawing velocity is in a simple 
ratio to that produced per unit rolling velocity (equation (16)). 
Thus 


Cl r~W Ci p 


.The derivative has been derived previously and is given 
in table II. 

Derivatives 6V r and C „ r . — When the wing yaws, the 
antisymmetrical pressure distribution which is indicated by 
equation (15) will produce unbalanced suction forces at tho 
right and left wing tips and thereby give rise to lateral forces 
and yawing moments. In addition, skin friction will con- 
tribute a yawing moment. 

It appears that a reasonable although approximate evalu- 
ation of the tip suction forces in yawing can bo obtained by' 
means of the correspondence of yawing with rolling as utilized 
previously in deriving equation (16). This procedure docs 
not satisfy the Kutta-Joukowski requirement in the sideslip 
component of the stream velocity in yawing; however, this 
theoretical deviation is likely to be very small in tho actual 
flow. On the basis of these considerations, tho induced suc- 
tion forces on the wing tips per unit yawing velocity will bo 
related in the ratio a/B 2 to the corresponding induced suction 
forces per unit rolling velocity which were derived previously 
(section _en titled “Derivatives Cy and The contri- 

butions of the tip suction forces to CV r and C %r are, therefore, 


and 


CV,— c Yp 


C — “ n 
°"ri ~B 2 


where Cy and C n are given in tablo II. 

The effect of skin friction on the yawing moment due to 
yawing is 

N t = cos aJ_ A J n 0 D 0 | [(F-ry^+r^-j^^r— ^jfi+y^dxdy 

where the first bracketed term expresses the squaro of the 
resultant local velocity and /S is the local angle of sideslip: 


0 



Eliminating second-order terms and terms corresponding to 
symmetrical drag forces and converting N t to coefficient 
form yields 


4<ri 3 




RESULTS AND DISCUSSION 


As noted in the preceding analysis, the nondimensional 
stability derivatives which are presented in table II were 
derived with reference to principal body axes with the origin 
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at point These results may be transformed by 

means of the equations in the last column, of table EL to 
apply to stability axes with the origin at an arbitrary dis- 
tance x et from the midchord point. The stability axes are 
shown in figures 2 (b) and are obtained by a rotation of the 
principal body axes (fig. 2 (a)) through an angle a; the origin 
is then shifted a distance x et along the new ar-axis. The 
conversion to stability axes was obtained by means of the 
transformation formulas given in reference 16, with the 
omission of relatively unimportant terms compared to unity, 
such as a*. 

The formulas for the derivatives given in table II with 
reference to principal body axes are shown plotted in figures 4 
and 5 against the parameter AB. (Derivative —C Xu and 
those derivatives equal to zero are omitted from the figures.) 
These curves show the variation of the stability derivatives 
with aspect ratio for constant Mach number. The varia- 
tion with Mach number for constant aspect ratio is not 



Fionsi 4.— Variation of supersonic lateral stability derivatives with aspect ratio-Mach 
number parameter. Derivatives with respect to principal body aies; thin flat rectangular 
wins: no dihedral. (See table n for conversion to stability ares.) 


directly indicated, although it can be determined from the 
curves. These data are shown in figure 4 for the lateral 
stability derivatives and in figure 5 for the longitudinal 
stability derivatives. The data in figures 4 and 5 are to be 
used in conjunction with the transformation formulas pre- 
sented in table II to evaluate the derivatives with respect to 
stability axes. In the evaluation of these derivatives, many 
of the terms are likely to be relatively small; therefore, the 
expressions will be noticeably simplified when such terms are 
neglected in the computations. 

The results of the present investigation have been derived 
on the assumptions of zero thickness and small disturbances. 
Potential flow is assumed except in the case of C Mr and C Xu 
in which skin friction is considered. The practical effects 
of the Kutta-JoukowsM requirements which are introduced 
at the wing tips in sideslip and yawing are not definitely 
known. On tbis basis, the data shown in figures 4 and 5 
are expected to apply in general to thin wing sections for 
small steady motions, motions with small accelerations, or 
oscillatory motions of low frequency in which boundaiy- 
layer effects are not expected to be important. The applica- 
bility of the present theory to Mach numbers in the vicinity 
of unity, very high Mach numbers, or for very low aspect 
ratios is uncertain. 

The data in figure 4 (a) show that at supersonic speeds the 
derivative WCija. is negative in contrast with the behavior 



(b) Yawing-moment-ooefficient derivatives. 
Fiatnti 4.— Continued. 
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AB 

(c) SIde-force-ooefflcient derivatives. 

Figuhe 4.— Concluded. 

at subsonic speeds where positive values would be obtained. 
This phenomenon was pointed out for the infinitely long 
wing in reference 7 and its physical significance elaborated 
upon. For stability axes, the formula for Ci r ' (table II) 
indicates that another reversal in sign to a positive value 
occurs as the Mach number is increased beyond approxi- 
mately 1.41 for typical rectangular wings. (Also see 
reference 7 for the infinitely long wing.) 

The suction force at the leading edge of rectangular wings 
vanishes at supersonic speeds. This factor should have an 
important influence on the derivatives C nf> and C„ r as super- 
sonic speeds are attained. In the case of C„J (stability 
axes), the results of the present analysis indicate that at 
supersonic speeds the sign of C n J will have positive values 
in many typical cases in contrast to negative values normally 
obtained at subsonic speeds. In the case of C„ r or C n ' , the 
loss of leading-edge suction tends to be compensated by 
the spanwise compressibility effects associated with super- 
sonic speeds. 

As noted previously in the analysis, the Kutta-Joukowski 
condition is unlikely to be satisfied along the trailing wing 
tip for a typically rounded wing tip at small angles of sideslip. 



Fiottbe 5. — Variation of supersonio longitudinal stability derivatives with aspoct railo-Macb 
number parameter. Derivatives with respect to principal body axes; thin flat rectangular 
wing. (See table II for conversion to stability axes.) 
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(d) c_„-c. I _+c'- v 


Therefore, the results for Ct p in figure 4 (a) are applica- 
ble where the Kutt a-Joukowski condition along the wing 
trailing edge is not satisfied. In order to determine the 
effect on C tp of satisfying the Kutt a-Joukowski condition 
along the trailing wing tip, the formula for C tg which meets 
this requirement was also obtained and is as follows: 

n ~ a ( 1 

\AB 3A 2 B* ) 

A comparison of this formula with the data for C l(J given 
in figure 4 (a) indicates that the effect of satisfying the Kutta- 
Joukowski condition along the trailing wing tip reduces 
negatively the values of C ^ from those obtained by neglect- 
ing the Kutta-Joukowski condition. For example for B= I 
and AB=4, when the Kutta-Joukowski condition along the 
trailing wing tip is neglected, C t .=— 0.083a; and when 
the Kutta-Joukowski condition is satisfied along the trailing 
wing tip, Ci fi = 0.146a. Thus, it is expected that when 
the sideslip angle becomes large, the dihedral effect — G lff 
should be reduced significantly because of the Kutta- 
JoukowsM condition along the trailing wing tip. 

The longitudinal stability derivatives in figure 5 refer to an 
axis whose origin is located at the midchord point. The data 
in figure 5 (c) for BC ma show that rectangular wings, with 
reference to this origin, have an increasingly unstable 
pitching moment with decreasing aspect ratio which cor- 
responds to a forward shift in the aerodynamic center. For 
infinite aspect ratio, the aerodynamic center is located at 
the midchord point or BC„ a = 0. If the aspect ratio is 
decreased to a value of 4 for a Mach number of 1.41, figure 5 
indicates a forward shift of the aerodynamic center of 
0.025 chord. With constant Mach number, the ratio 
BC m J—BC Za is obtained from figure 5 solely as a function of 
AB. These data indicate that with constant aspect ratio 
and increasing Mach number, the aerodynamic center will 
shift rearward. For an aspect ratio of 4, an increase in 
Mach number from 1.4 to 1.9 will shift the aerodynamic 
center rearward 0.01 chord. 

The derivative — C x& given in table II for infinite aspect 
ratio is negative which indicates negative damping or insta- 
bility. The ratio — C M J— C z& from table II gives the location 
of the center of pressure of the resultant lift contributed 
by a. By taking this ratio for infinite aspect ratio, the 

2 

center of pressure is found to be located at a point g c behind 

the leading edge. The negative damping produced by a, 
therefore, gives an unstable pitching moment for center-of- 

2 

gravity locations ahead of z=g c. These unstable tendencies 

caused by a are m inimiz ed by the effects of finite span and 

.M*4-l 

the instability due to C z - disappears entirely if AB ^ — ^ — . 


Fiocei 6.— Concluded. 
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TABLE I —DISTRIBUTIONS OP ADDITIONAL SURFACE VELOCITY POTENTIAL AND PRESSURE-DIFFERENCE COEFFICIENT CAUSED BY 

VARIOUS WING MOTIONS 


[Thin flat rectangular wing; no dihedral; AB^. 2; axes are shown in figs. 1 and 2] 
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TABLE II.— STABILITY DERIVATIVES OF THIN FLAT RECTANGULAR WINGS WITHOUT DIHEDRAL AT SUPERSONIC SPEEDS 
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Stability axes 


(^origin at point (^, 0, 0^ 

^origin at distance x cg measured positive ahead of midchord point x egi 0,0 jj 
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